Pseudo-random operators of the circular ensembles 
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We demonstrate quantum algorithms to implement pseudo-random operators that closely repro- 
duce statistical properties of random matrices from the three universal classes: unitary, symmetric, 
and symplectic. Modified versions of the algorithms are introduced for the less experimentally 
challenging quantum cellular automata. For implementing pseudo-random symplectic operators we 
provide gate sequences for the unitary part of the time-reversal operator. 
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The possibility of manipulating, transferring and stor- 
ing information in a way that preserves quantum coher- 
ence has led to a reexamination and extension of the pos- 
tulates of information processing ^ . This field of study, 
known as quantum information processing, can claim as 
triumphs the discovery of quantum algorithms that factor 
large numbers 0] , search databases quickly , and sim- 
ulate quantum systems efficiently ^4J. In addition, pow- 
erful communication and cryptographic protocols have 
been suggested based on the laws of quantum mechanics 

Q 

Random number generation is a basic component of 
classical information theory. Their quantum counter- 
parts, random quantum states and operators, likewise 
play a vital role in quantum information theory. Quan- 
tum communication protocols utilizing randomness in- 
clude saturation of the classical communication capacity 
of a noisy quantum channel by random states and su- 
erdense coding of quantum states via random operators 
Quantum computing protocols facilitated by random 
unitaries include quantum process tomography via a ran- 
dom operator fidelity decay experiment to identify types 
and strengths of noise generators 0- In addition, the 
amount of multi-partite entanglement in random states 
approaches the maximum at a rate exponential with the 
number of qubits in the system |8l|. 

Random quantum states, generated by applying a ran- 
dom operator to computational basis state, can also 
be used for unbiased sampling. When testing an algo- 
rithm, such as quantum teleportation, or a communica- 
tions scheme it is desirable to insure success for all pos- 
sible quantum states. This can be done via quantum 
rocess tomography, however this is extremely inefficient 
. Rather, one could test the likelihood of success with 
states drawn in an unbiased manner from the space of all 
quantum states, similar to sampling statistics in other 
contexts. 
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To capitalize on the above uses of random states and 
operators, it is necessary to efficiently implement ran- 
dom matrices on a quantum computer. This would ap- 
pear to be a daunting task considering that the num- 
ber of independent variables in a given random operator 
grows exponentially with matrix dimension. Neverthe- 
less, pseudo-random operators suggest that it may be 
possible to efficiently reproduce statistical properties of 
randomness on a quantum computer In this pa- 

per we extend the algorithm of Ref. [3 to produce 
pseudo-random operators from the universal random ma- 
trix classes with time-reversal symmetry. 

Pseudo-random operators from the universal classes 
with time-reversal symmetry may help a quantum com- 
puter simulate systems with time-reversal symmetries. 
These systems are especially important in the areas of 
quantum chaos 0| and decoherence , where the phys- 
ical system or environment to be studied tend to be mod- 
eled with time-reversal symmetry. We note that models 
of decoherence, specifically with random classical fields, 
have already been implemented on a nuclear magnetic 
resonance quantum information processor jl2l |. 

In addition, we extend our recent work and show 
that random operators from all three universal classes 
can be implemented on the less experimentally challeng- 
ing quantum cellular automata (QCA) architecture. This 
further demonstrates the usefulness of a QCA in the 
study of complex quantum evolution. For both architec- 
tures we show how to implement the unitary part of the 
time-reversal operator on a quantum computer which, 
for the QCA case, requires identifying a suitable, non- 
standard form of the time-reversal operator. 

Random matrices were first introduced by Wigner to 
describe the energy level spacings of large nuclei 
Since then, random matrices have functioned as a uni- 
versal model for a host of complex systems ranging from 
quantum dots to field theory X^. The circular ensembles 
of unitary matrices were introduced by Dyson [l^ as al- 
ternatives to the Gaussian ensembles of Hermitian matri- 
ces 0, ^3 • The three circular ensembles are the circular 
unitary ensemble (CUE) of arbitrary unitary matrices, 
appropriate for modeling systems without time rever- 
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sal symmetry, the circular orthogonal ensemble (COE) 
of symmetric unitary matrices, appropriate for systems 
having time-reversal invariance and integral spin or ro- 
tational symmetry, and the circular symplectic ensemble 
(CSE) of self-dual unitary quaternion matrices, appropri- 
ate for systems with time-reversal invariance, half-integer 
spin, and no rotational symmetries. Each universality 
class has properties unique unto itself. For example, the 
degree of level repulsion (the rate of change of nearest 
neighbor eigenangle spacings as the spacing goes to zero) 
is P(s) - s for COE, P{s) - for CUE, and P{s) - s'^ 
for CSE, where s is the nearest neighbor eigenangle spac- 
ing. Additionally, the distribution of eigenvector compo- 
nent amplitudes follow the xt distribution |l8l | 
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where y is the eigenvector component amplitude, and 
{y) is the mean value of y. The number of degrees of 
freedom, i/, is 1 for the orthogonal ensemble, 2 for the 
unitary ensemble, and 4 for the symplectic ensemble. 

The algorithm introduced in 7] to produce pseudo- 
random operators of arbitrary unitaries, CUE, consists of 
TO iterations of the n qubit gate: apply a random SU(2) 
rotation to each qubit, then evolve the system via all 
nearest neighbor couplings 7]. A random SU(2) rotation 
on qubit j of iteration i is defined as 18] 
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where the angles 0^ , and V'i are drawn uniformly from 
the intervals 



< < 27r < i^'y < 2tt, 



(3) 
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and 0^ = sin~^(^^ ) where is drawn uniformly from 
to 1. The nearest neighbor coupling operator at every 
iteration is 

Unnc = exp (^{7t/4) E ® '^l'"'^ : (4) 

where al is the z-direction Pauli spin operator. The ran- 
dom rotations are different for each qubit and each iter- 
ation, but the coupling constant is always 7r/4 to max- 
imize entanglement. After the m iterations, a final set 
of random rotations is applied. This algorithm has been 
shown, for up to 10 qubits, to implement operators with 
statistical properties extremely close to those expected of 
CUE with relatively few iterations [l^, Fig. [l] 

Random operators from the other two universal classes 
can be constructed from CUE operators. The goal of this 
paper is to demonstrate that these constructions allow for 
simple modifications of the pseudo-random algorithm de- 
scribed above to generate pseudo-random operators from 
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FIG. 1: Distributions of nearest neighbor eigenangle spac- 
ings, s (left), and eigenvector component amplitudes, y 
(right), for 100 realizations of 8-qubit, 60 iteration pseudo- 
random CUE operators compared to that expected for ran- 
dom unitaries of COE (dash), CUE (solid), and CSE (dash- 
dot). The nearest neighbor eigenangle spacing distribu- 
tion compares extremely well with the expected distribution 
Pcue{s) = (32s'^/7r^) exp(— 4s'^/7r) and the operators' eigen- 
vector component amplitude distribution almost exactly fol- 
lows Pcusiy) ~ exp(— J/), which is appropriate when (y) = 1 
and in the limit N —> oo. 



these classes. To draw a matrix, Ucoe from the space of 
all COE simply draw Ucue, from CUE and multiply it 
by its transpose 0, [H | 



Ucoe — Uqu^Ucue- 



(5) 



Pseudo-random generation of such an operator is read- 
ily done. First, implement the CUE matrix Ucue as 
above, retaining in memory (quantum or classical) the 
values of the 3n(m -I- 1) -I- 1 independent variables nec- 
essary to implement the operator, 3 for each rotation 
of n qubits for m + 1 rotations and 1 for the coupling 
strength. Next, implement the transpose of the operator 
by applying the transpose of each specific operation in 
reverse order of the original Ucu e generation. To imple- 
ment the transpose of the rotations apply the transpose 
of each individual qubit rotations. The transpose of the 
coupling operation is the same coupling C/,"f„j, = Unnc- 
In Fig. 121 we demonstrate that for 8 qubits and 60 it- 
erations the pseudo-random COE matrices generated in 
this way fulfill statistics of randomness by comparing the 
distributions P{s) and P{y) to that expected for COE. 

We now turn to the symplectic ensemble, representing 
systems with half-integer spin that are invariant under 
time-reversal. Following Mehta ^3 we define the anti- 
unitary time-reversal operator T = ZC , where C takes 
the complex conjugate and Z is unitary. The symplectic 
ensemble is characterized by an anti-symmetric Z, i.e. 
ZZ* = — 1 where * means non-Hermitian conjugate. We 
choose the representation such that Z is written Ii®!^® 
. . . ® ® Zn where Ij is the two-dimensional identity 
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FIG. 2: Distributions of nearest neighbor eigenangle spac- 
ings, s (left), and eigenvector component amplitudes, y 
(right), for 100 realizations of 8-qubit, 60 iteration pseudo- 
random COE operators. The nearest neighbor eigenangle dis- 
tribution compares very well with the distribution Pcoe{s) — 
(7rs/2) exp(— 7rs^/4) and the eigenvector component am- 
plitude distribution almost exactly follows PcoE{y) ~ 
(l/V2^)exp(-y/2). 

matrix and 

In this way, a symplectic unitary is defined by 

UcSE = -ZUcsE^ = UcSE- (7) 

and UcseZ is anti-symmetric unitary. 

As with the COE operators, drawing an operator from 
CSE can be done via CUE operators: draw Ucue and 
multiply by its time-reversal [l7ll2(J | 

UcsE = Uqu^Ucue, (8) 

where U^^jj^ = —ZU]jjj^Z . Using this construction 
pseudo-random CSE operators can be generated as fol- 
lows: run the pseudo-random operator algorithm to im- 
plement Ucue, apply Z via two rotations of the least 
significant qubit z = exp(— i(7r/2)cr2) exp(— i(7r/2)(T2:), 
where the ai are the Pauli matrices, Uqjj^ is imple- 
mented as explained in the COE case, apply Z. The 
negative sign is a global phase. 

Figure 121 shows the eigenangle and eigenvector element 
distributions for CSE pseudo-random operators. We note 
that matrices of the CSE exhibit Kramers' degeneracy 
so we digress to explain how the above distributions are 
determined. Kramers' degeneracy allows the following 
basis choice for CSE matrices ^3 

|l),T|l),|2),r|l),...|7V/2),T|7V/2), (9) 

where T is the time-reversal operator. The nearest neigh- 
bor eigenangle distribution uses only one of each degen- 
erate eigenangle. The eigenvectors corresponding to the 
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FIG. 3: Distributions of nearest neighbor eigenangles, s 
(left), and eigenvector component amplitudes, y (right), 
for 100 realizations of 8-qubit, 60 iteration pseudo-random 
CSE operators. The nearest neighbor eigenangle distribu- 
tion compares very well with the distribution Pcse{s) = 
(64/971)"^ exp(—64s^/97r) and the eigenvector component 
amplitude distribution almost exactly follows Pcssiy) = 
% exp(-2y). 

degenerate eigenangles can be written as 

|ei) - ci|l) + 5iT|l) + C2|2)+S2r|2)... 
r|ei) = -cJ|l) + cjr|l)-5;|2)+c*T|2)... (10) 

Any given diagonalization code will not necessarily out- 
put the above form for the two eigenvectors of a degen- 
erate eigenvalue, but superpositions of the two. Thus, as 
an invariant quantity to characterize the eigenvectors we 
use y = lcil^-|-|cip |lfl|. Using the above procedures, the 
distribution of nearest neighbor eigenangle and eigenvec- 
tor components for the generated pseudo-random CSE 
operators are those shown in Fig. |3| 

Classical cellular automata have been used to simu- 
late many complex classical systems from crystal growth 
to fluid flow |2l| . Thus, one may expect that quantum 
cellular automata (QCA) can be used to model complex 
quantum systems. Ref. [ij] demonstrates the implemen- 
tation of CUE pseudo-random operators. Here we extend 
that work to the other two random matrix classes. 

A QCA system is devised of k species of qubits in which 
all qubits of a species are addressed simultaneously and 
equivalently. Experimental flexibility is a primary moti- 
vation to explore implementations via QCA. Removing 
the need for localized external Hamiltonians can greatly 
ease hardware specifications for actual implementations 
of quantum information processing. A number of works 
have been devoted to exploring the universality of QCA 
architectures |3| but, despite the greater ex- 

perimental ease of QCA, relatively little work has been 
done to exploit the uniqueness of the QCA architecture 
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Previous work has shown that the pseudo-random al- 
gorithm appHed to a one-species QCA chain, such that 
all qubits undergo the same rotations, yields operators 
with eigenvalue and eigenvector distributions apprqpri- 
ate for CUE-type operators with mirror symmetries |l3| . 
The use of a two species QCA with alternating qubit 
species or the change of one nearest neighbor coupling 
constant (say from 7r/4 to 7r/5) is sufficient to break this 
symmetry. 

For a QCA COE pseudo-random operator, one ap- 
plies the pseudo-random operator algorithm, with all 
qubits of a species undergoing the same rotation, fol- 
lowed by its transpose, just as in the circuit architecture. 
The eigenvector component amplitude distribution for 8- 
qubit, one- and two-species COE operators is shown in 

fig. a 

To generate CSE pseudo-random operators requires 
the implementation of Z which above was done by in- 
dividually address the least-significant qubit. This oper- 
ation is illegal on a QCA system. Thus, we must find 
an appropriate non-standard representation of Z which 
allows all qubits of a species to be addressed equivalently. 

To find a representation of Z amenable to a QCA 
implementation, we recall that for a symplectic matrix, 
UcsE, the matrix A = UcseZ is antisymmetric unitary. 
For every antisymmetric unitary matrix there exists a 
unitary matrix W such that A = WZW'^ We de- 

fine our modified operator as Z' = VZV'^ . Since Z' Z'* = 
— 1, by definition of a symplectic matrix, {V*Y'V — 
V{V*Y' = ±1- Thus, V must be a symmetric or an- 
tisymmetric unitary. We can then define the antisym- 
metric unitary A' = UcseZ' = W Z'VV^ and, following 
Mehta [l^, we choose the unitary Ucue = {Z'W'y and 
generate symplectic matrices via Ucse = ~Z'Uqjj^Z'U . 
An example of a symmetric unitary operator, V ^ that al- 
lows the generation of CSE operators in the above fashion 
is the swap gate. This should not be surprising as the or- 
dering of qubits is completely arbitrary. An operator Z' 
that is appropriate for our purposes is the rotation z ap- 
plied to an odd number of qubits. Using this form of Z', 
pseudo-random symplectic operators (with mirror sym- 
metry if all coupling constants are equal) can be imple- 
mented on a one-specie QCA if there are an odd number 
of qubits. Similarly, two-species QCA implementations 
of pseudo-random symplectic operators can be achieved 
if one of the species consists an odd number of qubits. 
The eigenvector component amplitude distribution of 7- 
qubit one- and two-species QCA operators is shown in 
Fig. H 

In conclusion, we have demonstrated quantum algo- 
rithms for pseudo-random operators from the COE and 
CSE universal classes. As with the original CUE pseudo- 
random operators P], we provide evidence that suggests 
that these operators may be able to fulfill statistical 
properties of these ensembles with an efficient number 
of gates. Efficient performance of such operators could 
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FIG. 4: Distributions of eigenvector component amplitudes, 
y, for quantum cellular automata (QCA) based operators: 
(A) 100 realizations of 8-qubit one-specie COE, (B) 200 re- 
alizations of 7-qubit one-specie CSE, (C) 100 realizations of 
8-qubit two-species COE, (D) 200 realizations of 7-qubit two- 
species CSE. All distributions are for m = 40 iterations. 
The distributions of the two-species QCA operators are indis- 
tinguishable from those of random COE and CSE matrices. 
However, the one-specie QCA operator distributions deviate 
from the random distributions due to mirror symmetry. We 
note that this symmetry would be broken in any actual ex- 
perimental implementation due to unequal nearest-neighbor 
couplings. 



be useful in simulating various complex quantum sys- 
tems. Similar operators can also be implemented using 
the less experimentally demanding QCA. This is a fur- 
ther demonstration that a QCA system can be a useful 
tool in the study randomness. 
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